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Third  Year  Final  Report 

(1  July  1996  -  31  December  1997) 

Abstract 

We  have  concluded,  within  the  three  year  contract  (1  July  1994  -  30  June  1997, 
extended  until  31  December  1997),  the  investigation  of  the  dynamics  of  a  series  of 
elementary  chemical  reactions  of  atoms  (0(^P,*D)  and  N(^D))  and  radicals  (OH(^n)) 
relevant  to  the  local  atmosphere  of  the  Space  Shuttle  in  low  Earth  orbit.  The  technique 
employed  was  the  crossed  molecular  beams  scattering  method  with  mass  spectrometric 
detection.  Specifically,  during  the  third  year  of  the  research  project  we  have: 

1)  Extended  further  the  investigation  of  the  important  reaction  0(^D)+H2  by  performing 
measurements  of  reactive  differential  cross  sections  at  more  collision  energies, 
stimulated  by  recent  theoretical  developments  (i.e.  the  calculation  of  new  accurate 
potential  energy  surfaces  and  new  dynamical  predictions  on  these  surfaces). 

2)  Carried  out  the  investigation  of  the  important  reactions  of  electronically  excited  N 
atoms,  N(^D),  with  molecular  hydrogen,  acetylene  and  ethylene.  These  reactions  were 
studied  in  substitution  of  the  initially  planned  reactions  of  the  NH  radical. 

3)  Studied  the  reactions  0(^P,*D)+CH3l  to  explore  the  role  of  nonadiabatic  effects 
(inter-system  crossing)  in  the  dynamics  of  atomic  oxygen  reactions.  These  reactions 
were  studied  in  substitution  of  those  of  0(*D)+CH4  and  0(^D)+H20  which,  for 
technical  difficulties,  could  not  be  investigated. 

4)  Constructed  a  second  radio-frequency  discharge  beam  source  which  will  allow  us  to 
study  the  dynamics  of  radical-radical  reactions  in  the  next  future. 

During  the  overall  three  years  period  of  the  contract,  we  have  therefore  studied 
the  following  reactions  : 

1)  0H  +  C0->C02  +  H 

2)  OH  +  H2  ->  H2O  +  H 

3)  0('D)  +  H2  oh  +  H 

4)  0(^P,^D)  +  CH3l  ^I0  +  CH3 

5)  N(^D)  +  H2^NH  +  H 

6)  Nf D)  +  C2H2  HCCN  +  H 

7)  N(^D)  +  C2H4  CH3CN  +  H 

Reactions  (l)-(3),  and  presumably  also  reactions  (5)  and  (6),  occur  in  the  local 
atmosphere  of  the  orbiting  Space  Shuttle.  Their  reaction  dynamics  have  been  elucidated 
in  detail  by  measuring  reactive  differential  cross  sections  in  crossed  molecular  beam 
scattering  experiments.  It  has  been  shown  that  OH  radicals  fi-om  0(1d)+H2,  NH  radicals 
from  N(^D)+H2,  HCCN  radicals  from  N(2d)+C2H2,  and  H2O  molecules  from  OH+H2 
are  formed  highly  vibrationally  excited  and  therefore  give  rise  to  strong  emission  in  the 
infrared.  This  information  is  particularly  relevant  for  a  detailed  understanding  and 
modeling  of  the  local  atmosphere  of  the  orbiting  Space  Shuttle.  CO2  from  reaction  (1)  is 
instead  found  to  be  formed  translationally  very  excited,  rotationally  excited  and,  hence, 
vibrationally  cold.  The  study  of  reaction  (4)  has  shown  that  triplet-singlet  inter-system¬ 
crossing  occurs  with  very  high  probability  when  a  heavy  atom  such  as  iodine  is  involved. 
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In  most  cases  the  experimental  results  have  been  compared,  within  fruitful 
collaborations,  with  state-of-the-art  dynamical  calculations  (quasiclassical  and/or 
quantum)  on  ab  initio  potential  energy  surfaces. 

The  information  obtained  from  the  study  of  the  above  reactive  collisions  at  the 
microscopic  level,  besides  their  interest  from  a  fundamental  point  of  view,  will  enable  the 
world  scientific  community  and  the  U.S.  Air  Force  to  better  understand  the  local 
environment  of  measuring  platforms  in  space  and  may  also  be  useful  for  the  improvement 
of  the  SOCRATES  code,  aimed  at  assessing  the  effects  of  contamination  on 
measurements  aboard  spacecrafts  in  low  Earth  orbit. 

A  detailed  report  on  reaction  (1),  in  which  experimental  results  and  theoretical 
predictions  are  compared,  has  been  published  in  a  chapter  in  Advanced  Series  in 
Chemical  Physics,  vol.  6  “The  Chemical  Dynamics  and  Kinetics  of  Small  Radicals’’ 
(Publication  n.  1). 

A  detailed  combined  experimental  and  theoretical  study  of  reaction  (2)  has  been 
published  in  Chemical  Physics  (Publication  n.  2). 

A  comprehensive  publication  on  reactions  (3)  is  currently  in  press  in  the  Journal 
of  Chemical  Physics  (Publication  n.  3)  and  some  results  will  appear  as  Comments  to  the 
Faraday  Discussion  no.  108  on  “The  Dynamics  of  Electronically  Excited  States” 
(Publications  n.  4  and  5). 

The  development  and  characterization  by  Stern-Gerlach  magnetic  analysis  of 
intense  supersonic  beams  of  ground  and  electronically  excited  oxygen  and  nitrogen 
atoms,  0(^P,  ^D)  and  N('‘S,  ^D),  as  well  as  of  Cl(^P3/2,i/2)  has  been  published  in  Israel 
Journal  of  Chemistry  (special  issue  on  “Molecular  beams”)  (Publication  n.  6). 

A  publication  on  reactions  (4)  is  currently  in  preparation  (Publication  n.  7),  while 
a  first  partial  report  will  appear  as  a  Comment  to  the  Faraday  Discussion  no.  108  on 
“The  Dynamics  of  Electronically  Excited  States”  (Publication  n.  8). 

A  publication  is  in  preparation  on  reaction  (5)  in  collaboration  with  theoretical 
groups,  in  which  the  experimental  results  are  compared  with  dynamical  predictions  on  a 
high  quality  ab  initio  potential  energy  surface  (Publication  n.  9).  A  preliminary  report 
will  appear  as  a  Comment  to  the  Faraday  Discussion  no.  108  on  “The  Dynamics  of 
Electronically  Excited  States”  (Publication  n.  10). 

A  publication  is  also  in  preparation  on  reaction  (6)  (Publication  n.  11)  for  which, 
again,  a  preliminary  report  will  appear  as  a  Comment  to  the  Faraday  Discussion  no.  108 
on  “The  Dynamics  of  Electronically  Excited  States”  (Publication  n.  12). 

A  publication  is  being  planned  also  on  reaction  (7)  (Publication  n.  13). 

In  addition,  acknowledgments  to  AFOSR  and  EOARD  is  also  made  in  a  series  of 
publications  on  related  reactions  (Publications  14-19). 

Finally,  the  above  results  have  also  been  presented,  during  the  last  three  years,  in 
numerous  (more  than  ten)  Invited  Lectures,  as  well  as  in  many  poster  contributions,  at 
international  conferences. 
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reaction” 

Comment  to  Faraday  Discussion  no.  108  (University  of  Sussex,  Brighton,  U.K., 
Dec.  15-17,  1997),  in  press. 
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On  “Comparison  between  experimental  differential  cross  sections  and 
quasiclassical  trajectory  calculations  on  ground  and  excited  ab  initio  potential 
energy  surfaces  for  the  0(‘D)+H2  reaction” 

Comment  to  Faraday  Discussion  no.  108  (University  of  Sussex,  Brighton,  U.K., 

Dec.  15-17,  1997),  in  press. 

6.  M.  Alagia,  V.  Aquilanti,  D.  Ascenzi,  N.  Balucani,  D.  Cappelletti,  L.  Cartechini,  P. 
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Third  Year  Report 

(1  July  1996  -  31  December  1997) 

During  the  3^^  year  of  the  research  project,  the  plan  was  to  extend  the  same  type  of 
measurements  carried  during  the  l^t  and  2^^  year  on  the  reactions  0(1d)+H2^D2), 
OH+H2,  OH+CO,  and  N(2d)+C2H2  to  also  other  systems  of  relevance  in  the  local 
atmosphere  of  the  orbiting  space  shuttle.  Specifically,  we  have: 

1)  Completed  the  investigation  of  the  important  reaction  0(1d)+H2  by 
performing  measurements  of  reactive  differential  cross  sections  at  more  collision 
energies,  stimulated  by  recent  theoretical  developments  (i.e.,  the  calculation  of  a 
new  accurate  potential  energy  surface  (PES)  and  new  theoretical  predictions  on 
this  PES). 

2)  Studied  the  reactions  0(^P,  *D)  +  CH3I. 

3)  Characterized  in  detail  our  supersonic  beams  of  atomic  oxygen,  nitrogen  and 
chlorine  by  using  the  Stern-Gerlach  magnetic  analysis  technique. 

4)  Completed  the  investigation  of  the  important  reactions  of  electronically 
excited  N  atoms,  N(2d),  with  H2  and  acetylene. 

5)  Extended  the  study  of  N(^D)  reactions  to  N(^D)  +  ethylene. 

6)  Built  a  second  radio-fi-equency  (RF)  discharge  beam  source  to  study  the 
dynamics  of  radical-radical  reactions  in  the  next  fUture. 

The  reactions  0(^P,^D)  +  CH3I  were  studied  in  substitution  of  the  reactions 
0('D)+H20  and  0(^D)+CH4,  originally  planned,  but  which  could  not  be  investigated  due 
to  technical  difficulties.*  The  OfP,  *D)  +  CH3I  system  allowed  us  to  explore  the  effect  of 
electronic  excitation  and  nonadiabatic  effects  (inter-system  crossing)  in  the  dynamics  of 
atomic  oxygen  reactions,  a  subject  of  practical  as  well  as  fundamental  interest  (in  this 
respect  it  can  be  regarded  as  a  prototype  reaction). 

The  investigation  of  N(^D)  reactions,  originally  not  planned,  as  those  with  H2 
C2H2,  and  C2H4,  was  carried  out  to  compensate  for  the  fact  that,  due  to  technical 
difficulties  encountered  during  the  preliminary  experiments  aimed  at  the  generation  of 
supersonic  beams  ofNH  radicals  and  to  the  limited  remaining  amount  of  time,  we  could 
not  pursue  the  investigation  of  the  dynamics  of  the  reaction  NH  +  NO. 

The  results  obtained  during  the  3’^'*  year  of  the  contract  can  be  summarized  as 
follows. 


0(*D)  +  H2 

The  0(*D)  +  H2  ^  OH  +  H  reaction  has  attracted  over  the  years  a  large 
attention,  both  experimentally  and  theoretically,  due  to  its  fundamental  importance  in 
combustion  and  atmospheric  chemistry,  and  to  the  fact  that  it  constitutes  a  prototype  for 
reactions  involving  a  strongly  bound  intermediate.  For  some  time  this  reaction  has  been 
thought  to  be  a  typical  insertion  reaction  occurring  on  the  ground  state  PES 
(corresponding  to  the  ground  state  of  H2O).  But,  recently,  this  notion  has  been 
challenged  by  new  experimental  and  theoretical  evidence  from  our  and  other 
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laboratories.  Indeed,  during  the  last  two  years  an  unprecedented  renewed  interest  has 
arisen  around  this  system.  Experimentally,  it  has  come  under  investigation  at  Oxford^ 
where  by  using  polarized-laser-photolysis  and  polarized-sub-Doppler-LIF  detection  in  a 
bulb  state-selected  energy-averaged  differential  cross  sections  for  some  specific 
rovibrational  states  were  derived,  at  Taiwan^  where  by  using  Doppler  shifl;  techniques 
and  Doppler-selected  time-of-flight  methods  and  H(D)-atom  REMPI  detection 
differential  cross  sections  were  derived,  at  Boulder'*  where  new  rate  constants  were 
measured,  at  Heidelberg^  where  integral  reaction  cross  sections  were  determined,  and  at 
Berkeley®  where  differential  cross  sections  are  being  derived  by  using  ion-imaging 
techniques.  Theoretically,  the  reaction  has  been  investigated  at  Northwestern  University’, 
Princeton’,  Argonne’,  New  York*,  Chicago^,  Coimbra*®,  Birmingham**,  and  Madrid’  at 
the  level  of  both  electronic  structure  calculations  (to  derive  improved  potential  energy 
surfaces)  and  dynamical  calculations,  either  by  quasiclassical  trajectories  and 
approximate  quantum  methods.  A  central  issue  in  most  of  these  investigations  has  been 
the  role  and  contribution  of  excited  electronic  states.  This  point  has  been  in  particular 
addressed  by  Schatz  and  coworkers’’*’,  who  have  presented  detailed  studies  of  the 
0(*D)+H2  reaction.  However,  despite  the  vast  amount  of  work,  the  mechanism  and 
dynamics  of  the  reaction  are  still  not  well  understood. 

The  0(*D)  +  Ha  reaction  was  studied  at  two  collision  energies  (E0=3.O  and  4.0 
kcal/mol)  during  the  first  year  of  this  project.  Further  work  on  this  reaction  during  the 
second  year  was  prompted  by  (i)  the  observation  of  a  small  contamination  by  H2O  and 
OH  radicals  (and/or  by  and  natural  abundance)  in  the  0(1d)  beam  which  had 
to  be  accounted  for  and  (ii)  by  a  renewed  experimental  and  theoretical  interest  on  this 
fundamental  system.  The  amount  of  correction  for  0(^D)+H2  was  quantified  by 
measuring  the  elastic  scattering  at  m/e=16  from  H2  and  by  determining  the  ratio 
(m/e=17)/(m/e=16)  by  measuring  the  scattered  signal  at  m/e=17  and  m/e=16  from  He. 
The  correction  turned  out  to  be  significant  only  at  small  scattering  angles,  close  to  the 
0(^D)  beam;  the  maximum  value  being  of  about  10%  at  the  smallest  angle  (0=5°)  and 
then  rapidly  decreasing  to  zero  at  0r:15°.  Stimulated  by  very  recent,  accurate 
calculations  of  the  ground  and  excited  PES  and  by  dynamical  calculations  on  these 
PES’’*’,  we  have  carried  out  more  crossed  beam  experiments  at  different  collision 
energies  on  0(*D)+H2  and  on  the  isotopic  variant  0(*D)+D2  and  have  examined  the 
results  at  the  light  of  the  most  recent  theoretical  developments.  On  the  whole,  we  have 
measured  product  and  angular  velocity  distributions  in  the  laboratory  frame  for 
0(*D)+H2at  Ec=0.88,  1.9  and  3.0  kcal/mol  and  for  0(*D)+D2at  Ec=5.3  and  6.1  kcal/mol 
(see  Fig.  1),  and  have  derived  center-of-mass  (c.m.)  product  translational  energy  and 
angular  distributions  (The  reaction  with  D2  was  studied  using  a  beam  of  **0(*D),  which 
allowed  to  detect  the  OD  product  at  w/e=20,  which  is  free  of  any  contamination).*’  The 
c.m.  angular  distributions  are  nearly  backward-forward  symmetric,  but  exhibit  a  slight 
forward  bias  at  the  lowest  Ec  and  a  slight  backward  bias  which  increases  with  collision 
energy  at  the  other  Ec  (see  Fig.  2).  About  30%  of  the  total  available  energy  is  found  to 
be  channeled  into  product  translational  energy.  We  have  compared*’  the  experimental 
results  with  those  of  quasiclassical  trajectory  (QCT)  calculations  (in  collaboration  with 
the  group  of  Prof  J.J.  Sloan,  Univ.  of  Waterloo,  Canada)  on  the  ground  state  DIM 
(diatomic-in-molecule)  potential  energy  surface  of  Kuntz  et  al.*'*  These  calculations, 
carried  out  at  some  of  the  energies,  find  the  product  c.m.  angular  distributions  to  be 
backward-forward  symmetric  in  the  investigated  energy  range.  This  translate  into  lab 
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angular  distributions  which  exhibit  an  asymmetry  (because  of  the  c.m.-lab  Jacobian 
transformation)  opposite  to  that  experimentally  observed,  as  can  be  seen  from  Fig.  1. 
Interestingly,  very  recently  Schatz  and  coworkers*^  have  shown  that  the  DIM  ground 
state  PES  is  essentially  correct  and  have  shown  that  it  predicts  symmetric  angular 
distributions,  reflecting  an  insertion  dynamics,  similar  to  those  predicted  by  the  new 
RKHS  surface,  which  is  based  on  accurate  ab  initio  calculations.  Therefore,  we  attribute 
the  onset  of  a  preferential  backward  scattering  at  energies  >1.9  kcal/mol  to  the  onset  of  a 
second  micromechanism,  a  direct  abstraction  mechanism  evolving  on  one  or  both  of  the 
excited  (lA”  and  2 A’)  potential  energy  surfaces,  which  for  collinear  geometry  are 
degenerate  and  are  calculated  to  have  a  barrier  of  2.4  kcal/mol  (see  below).  Our  results 
are  in  line  with  recent  work  by  Liu  and  coworkers^  who  find  a  preferred  backward 
scattering  in  the  reaction  0(^D)+HD  at  Ec=4.55  kcal/mol  which  is  not  consistent  with  a 
simple  insertion  mechanism  on  the  ground  state  1  A’  PES  :  a  component  of  abstraction  in 
the  dynamics  was  invoked;  in  addition,  measurements  of  the  excitation  function^  have 
suggested  a  barrier  of  1.5  kcal/mol  for  the  abstraction  channel  in  the  reaction  0(^D)+H2. 
These  results  and  our  finding  of  a  preferred  backward  scattering  already  at  Ec=1.9 
kcal/mol  suggest  that  the  abstraction  barrier  in  the  new  ab  initio  PES  is  perhaps  slightly 
too  high. 

It  is  interesting  to  note  that  the  QCT  calculations  of  Schatz  et  al.^’*^  show  that  the 
ground  state  lA’  PES  gives  a  symmetric  angular  distribution,  while  the  excited  1  A”  PES 
gives  a  backward  peaked  angular  distribution  ;  the  global  (lA’  +  lA”)  angular 
distributions  for  0(‘D)+H2  at  Eo=2.7  and  5.0  kcal/mol  result  to  be  asymmetric  with  a 
backward  bias,  in  line  with  our  findings.  Also,  the  average  fraction  of  energy  in 
translation  is  calculated  by  Schatz  et  al.‘^  to  be  about  33%  at  Ec=5.0  kcal/mol,  which 
compares  well  with  the  value  of  about  30%  determined  in  our  experiments  at  slightly 
lower  Ec.  It  seems  that  the  new  PES  of  Schatz  et  al.  should  be  able  to  describe  the 
dynamic  of  the  0(^D)+H2  reaction  quite  well. 

Fig.  3  depicts  a  schematic  correlation  diagram  of  the  PES  which  shows  the  two 
possible  reaction  pathways  according  to  the  new  ab  initio  calculations’’*^:  on  the  ground 
state  PES  (*Z  and  lA*  in  Coov  and  C2v  geometry,  respectively)  insertion  of  0(’D)  into  the 
H2  bond  takes  place  without  energy  barrier  (for  linear  approach  the  barrier  is  0.7 
kcal/mol);  on  the  excited  lA”  PES  the  reaction  path  is  linear  with  a  barrier  of  2.4 
kcal/mol.  The  lA”  PES  correlates  adiabatically  with  ground  state  OH(’n)+H,  while  the 
other  Renner-Teller  component  of  *n,  the  2 A’  PES,  can  only  contribute  via  surface 
hopping  to  ground  state  products,  since  it  adiabatically  leads  to  OH(^E)+H.  From 
measurements  of  angular  and  velocity  distributions  of  the  OD  product  from  the  reaction 
0(*D)+D2  at  Ec=5.3  kcal/mol  we  have  derived  the  c.m.  angular  distribution  depicted  in 
Fig.  4.  This  is  compared  with  the  results  of  QCT  calculations  by  Schatz  et  al.  on  both  the 
ground  state  lA’  and  excited  state  lA”  surfaces. As  can  be  seen,  the  lA’  PES  gives  an 
angular  distribution  which  is  backward-forward  symmetric  and  reflects  an  insertion 
dynamics,  while  the  lA”  PES  gives  a  strongly  backward  peaked  angular  distribution  as 
one  would  expect  for  a  direct  abstraction  process.  Assuming  that  the  reaction  proceeds 
adiabatically  on  both  surfaces,  a  summed  contribution  (shown  in  Fig.  4  as  lA’+lA”)  is 
obtained  which  exhibits  a  backward  biased  angular  distribution  in  fairly  good  agreement 
with  experiment.  The  role  of  the  other  excited  2 A’  PES  (see  Fig.  3)  needs  still  to  be 
taken  into  account,  but  it  is  expected  to  contribute  considerably  less  than  lA”  and  its 
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angular  distribution  is  likely  to  be  relatively  flat  versus  0.'^  It  appears  that  dynamical 
calculations  which  take  into  account  both  the  ground  state  surface  (which  has  no  energy 
barriers  to  insertion)  and  the  excited  lA”  surface  (which  has  a  barrier  of  2.4  kcal/mol  to 
collinear  abstraction)  are  able  to  account  nearly  quantitatively  for  the  experimental 
backward  peaked  OD  angular  distribution  measured  for  the  0(*D)+D2  reaction  at  Ec=5.3 
kcal/mol.  The  abstraction  pathway  is  shown’’*^  to  give  highly  vibrationally  excited  OH 
product,  which  summed  to  the  vibrational  distribution  obtained  from  the  insertion 
pathway  leads  to  theoretical  results  which  are  in  line  with  experimental  observation  from 
infrared  chemiluminescence  experiments. 

Although  exact  quantum  scattering  calculations  need  to  be  performed  for  an 
accurate  assessment  of  the  new  PES,  we  feel  that  the  most  recent  theoretical  work  has 
led  to  an  improved,  though  still  not  fully  quantitative,  comparison  with  experiment  and 
then  we  have  now  a  better  understanding  of  the  dynamics  of  the  0(*D)3-H2  reaction. 

We  can  conclude  that  the  0(^D)+H2  reaction,  which  occurs  in  the  local 
atmosphere  of  the  orbiting  space  shuttle  and  is  also  the  main  source  of  OH  radicals  in  the 
stratosphere,  takes  place  through  a  direct-insertion  mechanism  at  low  collision  energies 
(<2  kcal/mol),  forming  a  highly  excited  H2O  molecule  intermediate,  which  is  very  short 
lived  and  dissociates  within  a  few  bending  vibrations,  before  the  energy  can  randomize  ; 
the  OH  angular  distribution  is  symmetric.  At  energies  above  2  kcal/mol  a  direct 
abstraction  mechanism  on  the  excited  PES  starts  to  contribute  producing  OH  highly 
vibrationally  excited  (v’=3  and  4  are  predicted’**^  to  be  the  most  populated  states  of  OH) 
and  leading  to  a  preferred  backward  scattering  of  the  OH  product  (with  respect  to  the 
direction  of  the  incoming  0(^D)).  The  high  degree  of  vibrational  excitation  of  the  OH 
product,  up  to  the  limit  of  energy  conservation,  clearly  indicates  that  the  OH  radical 
formed  from  the  0(^D)+H2  reaction  is  optically  very  active  and  therefore  strong 
emission  in  the  infrared  is  to  be  expected.  This  needs  to  be  taken  into  account  when 
modeling  the  local  atmosphere  of  the  Space  Shuttle  orbiting  around  the  Earth. 


0(''P,  ’D)  +  CH3I 

The  reactions  of  atomic  oxygen,  both  in  the  ground-state  0(^P)  and  the  first 
electronically  excited  state  0(^D),  with  halogenated  compounds  are  of  interest  in 
determining  the  impact  of  anthropogenic  surface  release  of  halogen  containing  molecules 
on  the  atmosphere  and  especially  on  the  ozone  natural  balance.  They  are  also  of 
relevance  in  the  combustion  chemistry  of  halogenated  compounds,  as  for  instance  in  their 
use  as  fire  extinguishers.  While  the  active  forms  of  chlorine  and  bromine  contribute  to 
stratospheric  ozone  depletion,  it  has  been  suggested  that  the  tropospheric  ozone  balance 
may  be  affected  significantly  by  iodine  compounds.  Recently,  it  has  been  proposed  that 
iodine  could  account  also  for  the  observed  depletion  of  lower  stratospheric  ozone  (below 
20  km,  altitudes  where  chlorine  and  bromine  are  not  very  effective  for  ozone  destruction 
in  contrast  to  iodine)  at  midlatitudes. The  main  source  of  atmospheric  iodine  is  natural 
(methyl  iodide,  chloroiodo-methane  and  diiodomethane  are  massively  produced  as 
metabolic  by-products  of  marine  biota),  but  now  iodine  containing  molecules  are  also 
considered  as  potential  substitutes  for  halon  fire  suppressants.  Although  CH3I  has  a 
short  lifetime  in  the  atmospheric  environment  due  to  its  rapid  photo-decomposition,  it  is 
interesting  to  probe  its  reactivity  with  atmospheric  atomic  and  radical  species,  such  as 
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0(^P)  and  also  0(^0)  which  is  abundant  and  very  active  in  the  ozonosphere.  Moreover, 
the  reaction  can  be  regarded  as  prototypical  of  the  reactions  of  0(3p,^D)  with  simple 
iodine  containing  molecules  and  their  study  allows  to  investigate  the  role  of  electronic 
excitation  in  reaction  dynamics  and  to  explore  the  possible  role  of  inter-system  crossing 
(ISC)  between  triplet  and  singlet  potential  energy  surfaces  (see  Fig.  5). 

Thermodinamically  allowed  products  for  the  reaction  0(^P)  +  CH3I  are: 


0(3p)  +  CH3I 


^  IO  +  CH3  (la) 

->  OH  +  CH2I  (lb) 

H  +  I  +  CH20(1c) 
^  I  +  CH3O  (Id) 

^  HI  +  CH2O  (le) 


while  the  list  is  even  larger  for  0(1d)+CH3I,  including  also  formation  of  HOI+CH2  and 
ICH2O+H. 

Measurements  of  the  absolute  rate  constants  between  213  and  364  K  and  product 
yields  at  298  K  have  been  recently  reported  for  the  0(^P)+CH3l  reaction,'*  The 
IO+CH3  and  OH+CH2I  channels  were  found  to  account  for  44%  and  1 6%  of  the  global 
rate  constant,  respectively;  product  yields  for  the  minor  channels  are  reported  to  be  ~7% 
for  H  formation,  <3%  for  CH3O  formation,  and  <5%  for  HI  formation.  While  in  the 
similar  reaction  0(^P)+CF3l  the  formation  of  10  is  largely  the  dominant  reaction 
channel,  the  observation  of  numerous  products  suggest  that  the  0(^P)+CH3l  reaction  is 
much  more  complex  and  occurs  on  multiple  potential  energy  surfaces.  No  dynamical 
information  was  available  on  this  system.  The  rate  constant  for  the  reaction  0(^D)+CH3l 
has  never  been  determined;  however,  recent  kinetic  studies  provided  room  temperature 
rate  constants  for  the  analogous  reactions  0(^D)+CH3C1  and  0(^D)+CH3Br  which  are 

1.9-10“^®  and  1.8- 10“^*^  cm^  molecule”^  s“^,  respectively.'^  The  10  forming  channel  is 
expected  to  be  the  dominant  reaction  pathway  also  for  the  0('D)  reaction, 

A  very  interesting  issue  is  the  extent  of  the  coupling  between  the  triplet  and 
singlet  potential  energy  surfaces  over  which  the  reaction  can  evolve.  A  crossing  between 
triplet  and  singlet  PES  always  occurs  in  0(*P,'D)  reactive  systems  since  0('D)  is  higher 
in  energy  than  0(*P)  by  about  2  eV  and  the  most  stable  intermediate  is  always  a  singlet 
species  (see  Fig.  5). 

Critical  for  the  study  of  the  0(^P,^D)+CH3l  reactions  by  the  crossed  molecular 
beam  scattering  method  has  been  our  capability  of  generating  intense  and  continuous 
supersonic  beams  of  oxygen  atoms  containing  0(3p)  and  a  small  percentage  of  0(^0).^“ 
We  have  carried  out  experiments  at  two  different  collision  energies  (E^),  13.4  and  15.5 
kcal/mol.  Under  our  experimental  conditions  we  have  been  able  to  observe  only  the  10 
product  (the  major  reaction  pathway'*)  from  the  reactions: 

0(*P)  +  CH3I 10  +  CH3  AH°o=-0.7±2.0  kcal/mol  (a) 

0('D)  +  CH3I  ->  10  +  CH3  AH°o=  -46,3±2.0  kcal/mol  (b) 


The  natural  abundance  of  in  the  CH3I  beam  produces  a  strong  interference  at  the 
mass  of  the  10  product,  making  the  detection  of  10  difficult  due  to  elastically  scattered 
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l^CH3l  contamination;  therefore,  to  overcome  this  complication  (which  is  severe  at 
angles  close  to  the  CH3I  beam),  we  used  isotopically  labeled  ^^02  as  precursor  of 
atomic  oxygen  and  detected  the  10  product  at  m/e=145.  Product  detection  for  reaction 
channels  (lb)  and  (Ic)  are  strongly  kinematically  unfavored,  while  product  detection  for 
channels  (Id)  and  (le)  are  hindered  by  the  dissociative  ionization  in  the  ionizer  of 
elastically  scattered  CH3I. 

From  10  product  angular  and  velocity  distributions^'  (see  Figs.  6-8),  the 
dynamics  of  the  reactions  (a)  and  (b)  were  characterized.  While  10  coming  from  the 
0(*D)  reaction  arises  from  the  decomposition  of  a  long-lived  complex  (the  dynamics  is 
similar  to  that  recently  found^^  for  the  related  reaction  0(*D)+HI->I0+H),  10  from  the 
0(^P)  reaction  is  found  to  be  formed  by  two  competing  micro-mechanisms  :  direct 
reaction  over  the  triplet  PES  giving  backward  scattered  10,  and  ISC  from  the  triplet  to 
the  singlet  PES  with  formation  of  a  bound  CH3IO  intermediate,  which  decomposes 
similarly  to  the  0('D)  adduct,  giving  a  backward-forward  symmetric  10  angular 
distribution  in  the  c.m.  frame.  The  best-fit  c.m.  angular  and  translational  energy 
distributions  at  Ec=13.2  kcal/mol  are  shown  in  Fig.  9.  In  Figs.  6-8  the  singlet  and  the  two 
triplet  contributions  to  the  best-fit  lab  angular  distributions  are  shown. 

The  resulting  weight  of  the  0(^D)  contribution  to  the  total  cm.  function  is  -35% 
of  that  of  0(3 P);  this  finding  is  not  surprising,  since  the  oxygen  beam  largely  consists  of 
0(3p),  but  the  rate  constant  of  0(^0),  which  is  not  known  but  can  be  estimated  from  the 
two  values  relative  to  the  analogous  reactions  0(1d)+CH3C1  and  0(^D)+CH3Br, 
should  be  about  one  order  of  magnitude  larger. 

The  fractions  of  reactive  encounters  which  undergo  ISC  is  found  to  be  about 
98%  of  the  total  0(^P)  reactive  collisions.  ISC  has  been  suggested  by  Grice  et  al.^^  to 
occur,  but  in  smaller  percentage,  in  the  related  0('*P)+RI  reactions  (R=C2H5,(CH3)2CH, 
(CH3)3C,  CF3CH2,  and  CF3)  and  has  been  correlated  with  the  effect  of  charge  transfer 
interaction  on  Renner  Teller  splitting  of  the  triplet  energy  surface.  Recent  ah  initio 
calculations  on  the  singlet  and  triplet  PES  of  the  0(^P,'D)+CH3I  system  support  the  role 
of  ISC  for  bent  geometry.^''  Our  results  indicate  that  ISC  is  the  dominant  pathway,  while 
the  abstraction  mechanism  accounts  for  only  a  few  percent  of  the  total  0(^P)  reaction  to 
form  10.  The  fact  that  the  latter  micro-mechanism  is  characterized  by  a  sharply  backward 
peaked  angular  distribution  (see  Figs.  6,  7  and  9)  suggests  that  only  for  strictly  collinear 
O-I-CH3  geometry,  within  a  narrow  cone  of  acceptance,  the  reaction  can  proceed 
directly  by  abstraction  on  the  triplet  PES.  The  direct  abstraction  contribution  is  found  to 
increase  slightly  as  the  collision  energy  is  raised  to  15.5  kcal/mol.  In  this  respect,  the 
dynamics  of  0(3p)+CH3l— >I0+CH3  appears  to  be  very  different  from  the  dynamics  of 
the  related  reactions  0(3p)-f-CF3l->I0+CF3  and  F(2p)-i-CH3l-^IF+CH3,  where  only  the 
complex  formation  mechanism  was  observed  to  occur.^^  It  would  be  interesting  to  have 
reliable  triplet  and  singlet  PES  for  the  system  and  to  carry  out  QCT  calculations  with 
surface  hopping  to  simulate  the  experimental  results  and  learn  more  about  the  role  of 
electronically  non  adiabatic  effects  in  oxygen  atom  reactions. 


Magnetic  analysis  of  supersonic  beams  of  O,  N  and  Cl  atoms 
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We  have  produced  supersonic  beams  of  O,  N,  and  Cl  atoms  and  of  metastable  N2 
and  O2  molecules  in  our  high-pressure  radio-frequency  discharge  beam  source  and  have 
characterized  them  by  coupling  velocity  selection  with  magnetic  analysis  in  the 
transmission  mode.^®  The  novelty  of  generating  supersonic  beams  of  N  atoms  of 
sufficient  intensity  to  carry  out  crossed  beam  reactive  scattering  experiments  was 
underlined  in  the  2"'*  year  report  and  is  exemplified  in  this  report.  The  results,  which  are 
described  in  detail  in  ref  26,  can  be  summarized  as  follows. 

We  have  found  that  oxygen  atoms  are  produced  mainly  in  the  electronic  ground 
state  ^Pj  and  a  substantial  relaxation  within  the  fine  level  structure  occurs  during  the 
supersonic  expansion.  An  upper  limit  of  10%  is  derived  for  the  relative  concentration  of 
atoms  in  the  first  electronically  excited  *D  state  in  the  beam. 

In  the  case  of  N  atoms,  not  only  the  ground  electronic  state  ‘‘S,  but  also  the 
metastable  and  states  are  formed  in  the  discharge  plasma ;  in  addition,  the 
relaxation  processes  are  ineffective  and  atoms  in  the  beams  have  a  high  concentration  of 
metastable  states  ('‘S  :  72%,  :  21%,  ^P  :  7%)  (see  Fig.  10).  Since  0(1D)  and  N(2D) 

atoms  are  metastable,  they  will  survive  up  to  the  interaction  region,  and  this  allows  us  to 
study  their  reaction  dynamics  in  crossed  molecular  beam  experiments. 

Cl  atoms  are  produced  exclusively  in  the  ground  electronic  ^Pj  state  (J=3/2,l/2), 
since  excited  electronic  states  are  energetically  inaccessible  at  the  plasma  temperatures. 
Atoms  in  the  excited  fine  level  ^Pi/2  are  produced  in  high  concentration  (about  15-20%) 
in  the  discharge  and  scarcely  relax  during  the  supersonic  expansion. 

We  have  shown  that  an  estimate  of  the  relative  population  of  electronic  states  in 
the  beam  on  the  basis  of  plasma  temperatures  or  final  translational  temperatures  can  led 
to  incorrect  conclusions. 

Finally,  magnetic  analysis  has  also  been  used  to  estimate  the  concentration  of 
molecular  oxygen  in  the  first  electronically  excited  a*Ag  state,  which  lies  «1  eV  above  the 
ground  state  and  has  a  lifetime  of  »44  minutes  and  can  be  present  in  the  beam  either 
because  of  incomplete  dissociation  or  because  of  atomic  recombination.  A  lower  and 
upper  limit  of  15%  and  40%,  respectively,  for  the  concentration  of  02(a^Ag)  in  the  beam 
was  derived  (see  Fig.  1 1),  Metastable  N2(A  and  N2(A’  which  have  an  energy 
of  6.3  eV  and  ^9.8  eV  respectively  over  the  ground  state,  have  also  been  found  to  be 
present  in  the  N  beams. 


NH  beam  generation 

Test  experiments  have  been  carried  out  to  ascertain  the  most  efficient  way  to 
generate  a  supersonic  beam  of  NH  radicals.  Unfortunately,  RP  discharge  in  dilute  (about 
2.5%)  mixtures  of  NH3  in  He  did  not  lead  to  NH  beams  of  sufficient  intensity  to  allow 
the  study  of  the  reaction  : 

NH  +  NO  ->  NNO  +  H  AH°^,=  -  34.7  kcal/mol 

The  goal  was  to  generate  a  beam  of  ^  ^NH  to  study  the  isotopic  variant  reaction 
i5nh  +  no  ^  *^NNO  +  H 
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by  detecting  N2O  at  a  mass-to-charge  ratio  m/e=45,  for  which  our  mass  spectrometric 
detector  has  an  inherent  background  of  a  factor  about  20  lower  than  for  m/e=AA,  but  due 
to  shortage  of  '^N2  extensive  tests  were  precluded.  Of  course,  these  attempts  will 
continue  in  the  next  future.  In  particular,  we  think  that  a  discharge  in  dilute  N2/H2 
mixtures  in  He  (here,  NH  would  originate  from  the  reaction  N(2d)+H2)  may  be  more 
efficient,  but,  again,  shortage  of  ^^N2  precluded  the  pursuit  of  this  investigation  for  the 
moment.  The  effort  was  therefore  concentrated  towards  the  completion  of  the  studies  of 
N(^D)  reactions  already  undertaken. 


N(2d)  +  H2 

The  reaction  of  ground  state  atomic  nitrogen  N(4s)  with  molecular  hydrogen  is 
considerably  endoergic  (AH=21.7  kcal/mol),  while  that  of  electronically  excited  N(2d)  is 
strongly  exoergic: 

N(2d)  +  H2  ->  NH(v')  +  H  AH°o=  -  33.2  kcal/mol 

and  may  lead  to  vibrationally  excited  imidogen  radicals,  in  analogy  to  the  corresponding 
0('d)  +  H2  — >  OH(v')  +  H  reaction.  The  room  temperature  rate  constant  for  the 
reaction  is  reported  to  be  about  5xl0"^2  Qn^3  rnolec'^s"^.  (the  reaction  of  N(^P)  is 
known  to  be  about  two  order  of  magnitude  slower^’  than  that  of  N(^D)). 

We  had  investigated  the  (technically  easier)  isotopic  variant: 

'^N(2d)  +  D2  ^  ^'ND(v')  +  D 

at  a  collision  energy  of  5.1  kcal/mol  (see  2"'’  year  report).  At  this  energy  the  ND  radical 
can  be  populated  up  to  v'=6  on  energetic  grounds.  The  angular  and  TOP  data  could  be 
fitted  with  a  symmetric  c.m.  angular  distribution  and  using  a  translational  energy 
distribution  reflecting  about  30%  of  the  available  energy  released  as  product  recoil 
energy.  This  indicates  that  the  ND  product  should  be  highly  internally  excited.  It  was 
inferred  that  the  reaction  proceeds  through  insertion  of  N(2d)  into  the  D-D  bond  to  form 
a  very  short-lived  NH2  radical  (amidogen)  which  dissociates  rapidly  after  a  few  bending 
vibration,  analogously  to  the  formation  of  highly  excited  OD  coming  from  the 
0(1d)+D2  reaction  through  a  short  lived  H2O  intermediate. 

During  the  past  year,  we  have  carried  out  new  angular  and  velocity  distribution 
measurements  at  a  lower  collision  energy  of  £^=3.8  kcal/mol  to  complement  the  data 
previously  taken  at  £^=5.1  kcal/mol  (see  Figs.  12  and  13).  It  is  found  that  also  at  this 
lower  £c  the  NH  product  is  formed  strongly  internally  excited  (about  30%  of  the  total 
available  energy  is  found  to  be  channeled  in  translation),  consistently  with  the  results  of 
recent  spectroscopic  investigation  of  this  reaction  by  Umemoto  and  coworkers. The 
c.m.  angular  distributions  at  £c=3.8  kcal/mol  (see  Fig.  14)  and  5.1  kcal/mol  are 
backward-forward  symmetric  and  reflect  an  insertion  dynamics  presumably  via  the 
ground  state  P£S  correlating  with  the  ground  state  ffBi  of  NH2.  The  absence  of  a 
preferential  backward  scattering  with  respect  to  the  similar  reaction  0('D)+D2  at  £c=5.3 
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kcal/mol*^  suggests  that  a  direct  abstraction  mechanism  does  not  occur,  at  least  in  this 
energy  range. 

Theoretical  calculations  of  the  experimental  observables  are  being  performed 
using  the  quasi  classical  trajectory  method  by  Schatz  and  coworkers^^  on  a  new  ab  initio 
PES  calculated  by  Harding,  Schatz  and  Rabitz.^**  Fig.  15  shows  the  schematic  energy 
level  and  correlation  diagram  for  the  reaction  N(^D)  +  H2  — >NH  +  H  according  to  the 
new  accurate  ab  initio  calculations  of  the  PESs.^**  The  insertion  pathway  on  the  ground 
state  surface  leading  to  the  l^Bi  state  of  NH2  is  found^®  to  have  a  barrier  of  about  1.7 
kcal/mol  in  C2V  geometry  and  of  about  4.7  kcal/mol  for  collinear  approach  (Coov 
geometry).  These  values  can  be  contrasted  with  corresponding  values  of  0  and  0.7 
kcal/mol,  respectively,  for  the  ground  state  PES  of  0('D)+H2.’’’^  The  first  two  excited 
PES  ^Ai  and  ^62  do  not  correlate  with  ground  state  NH(X^E'),  but  with  NH(a'A)  and 
NH(b^Z‘"),  respectively,  so  neither  can  contribute  adiabatically  to  the  formation  of  ground 
state  NH+H.  For  collinear  approach  the  first  excited  PES  has  a  large  barrier  (>  20 
kcal/mol)  and  does  not  correlate  with  ground  state  NH,  which  is  quite  different  to  what 
happens  in  the  related  0(‘D)+H2  reaction,  where  the  first  excited  PES  'FI  (in  collinear 
geometry)  has  a  barrier  of  2.4  kcal/mol  and  correlates  with  ground  state  OH.  These 
theoretical  calculations  suggest  that  only  insertion  plays  a  role  in  N(^D)+H2  at  thermal 
energies,  in  contrast  to  0(‘D)+H2.  The  0(‘D)+H2  and  N(^D)+H2  systems,  indeed, 
although  share  a  lot  of  similarities  (see  Fig.  16),  do  exhibit  considerable  differences  in  the 
topology  of  the  PES  (see  Fig.  3  and  15).  The  absence  of  a  backward  scattered 
contribution  in  the  angular  distribution  of  N(2d)+D2  with  respect  to  0(^D)+D2  can  be 
even  appreciated  directly  in  the  laboratory  angular  distribution  data  (see  Fig.  17) :  a  clear 
indication  that  abstraction  does  not  play  a  role  in  N(^D)+D2.  The  study  of  the  various 
surfaces  involved  in  this  reaction,  currently  in  progress  in  Schatz’ s  group, is  showing  a 
rather  complex  coupling  scheme  between  the  ground  state  PES  and  the  first  two  excited 
PES,  which  could  lead  to  a  significant  contribution  to  reaction  via  the  excited  state 
surfaces  even  at  fairly  low  collision  energies.  In  fact,  insertion  can  come  not  only  from 
the  ground  state  PES,  but  also  from  the  ^Ai  PES  which  is  adiabatically  connected  to  the 
^62  state  via  non-C2v  geometry  (a  conical  intersection  exists  between  ^Ai  and 
surfaces)  and  the  ^62  surface  has  a  modest  barrier  of  only  2.7  kcal/mol  (for  C2V 
geometry)  followed  by  a  fairly  deep  well  (see  Fig.  15). 

QCT  calculations  were  performed  by  Schatz  et  al.^^’^'  at  the  collision  energies  of 
the  experiments  on  the  ground  state  PES  and  calculated  product  angular  and 
translational  energy  distribution  were  compared  with  the  experimentally  derived 
quantities.  Quite  good  agreement  is  found  for  the  angular  distributions  as  can  be  seen  in 
Fig.  14  for  N(^D)+D2  at  Ec=3.8  kcal/mol  and  very  good  accord  exists  for  the 
translational  energy  distributions.^^’ Under  our  experimental  conditions  ND  can  only  be 
formed  in  the  ground  state,  but  raising  the  collision  energy  sufficiently  also  allows  for 
formation  of  ND(a*A)  (see  Fig.  15):  this  will  be  explored  in  future  experiments  together 
with  calculations  of  the  excited  ^Ai  PES  correlating  with  ND(a^A)  which  are  currently  in 
progress. 

Recently,  the  dynamics  of  the  N(^D)+H2  reaction  has  also  been  studied  by  QCT 
methods  on  a  previous  ab  initio  PES^^  and  the  results  are  quite  different  (backward 
scattered)  since  the  abstraction  barrier  on  this  surface  is  lower  than  the  insertion  barrier. 
Very  recently,^^  this  surface  has  been  modified  to  make  the  insertion  barrier  lower,  and 
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the  new  results  are  very  similar  to  what  obtained  on  the  new  PES.^'  Dynamical 
calculations  on  this  system  have  been  carried  out  by  quantum  scattering  methods  (QM) 
for  total  angular  momentum  J=0  using  several  of  the  available  potential  energy  surfaces, 
including  the  new  one  discussed  here.  The  results  are  generally  quite  similar  to  the  QCT 
results,  indicating  that  the  QCT  results  should  be  accurate. 

A  common  feature  to  the  N(^D)+H2  and  0(’D)+D2  reactions  is  that  the  radical 
product  is  highly  vibrationally  excited,  that  is,  infrared  active.  Therefore,  both  OH  and 
NH  emission  in  the  infrared  will  occur  in  an  environment  where  0(^D)  and  N(2d)  are 
formed  in  the  presence  of  H2,  as  it  can  occur  in  the  local  atmosphere  of  the  orbiting 
Space  Shuttle. 


N(^D)  +  C2H2 

In  the  2"'*  year  report  the  results  of  our  crossed  beam  investigation  of  the  reaction 

N(2d)  +  C2H2  HCCN  +  H  AH°o-  -  55  kcal/mol 

were  described.  Formation  of  the  HCCN  radical  from  this  reaction  has  been  suggested  to 
be  the  first  step  in  the  formation  of  cyanogen,  dicyanogen,  HCN,  etc.  in  the  atmosphere 
of  Titan. The  HCCN  radical  has  also  been  observed  in  interstellar  space^^  and  has 
attracted  much  experimental  and  theoretical  attention  over  the  years. Fig.  18  depicts  a 
schematic  energy  level  and  correlation  diagram  for  the  N(''S,  ^D,  ^P)  +  C2H2  system. 

From  product  angular  and  velocity  distributions  at  m/e=38  and  39,  which  were 
found  to  be  identical,  we  concluded  that  HCCN  formation  was  the  observed  reaction 
product,  while  the  slightly  more  exoergic  channel  leading  to  CCN  +  H2  does  not  occur 
appreciably.  The  HCCN  lab  angular  and  time-of-flight  distributions  at  £^=9.5  kcal/mol 
are  shown  in  Fig.  19  and  20,  respectively.  The  continuous  lines  represent  the  best-fit 
using  a  separable  form  of  the  center-of-mass  frame  product  flux  distribution,  m 
,E)=T(fi)P(E).  The  c.m.  angular  distribution  was  found  to  be  symmetric  at  Ec=3.1 
kcal/mol  and  somewhat  forward  peaked  at  the  higher  Ec  of  9.5  kcal/mol.  From  this  it  was 
deduced  that  the  reaction  proceeds  through  a  long-lived  complex  intermediate 
(presumably  the  cyanomethyl  radical  HgCCN)  at  low  E^  and  through  an  osculating 
complex  (i.e.  a  complex  having  a  lifetime  comparable  to  its  rotational  period)  at  the 
higher  The  cyanomethyl  radical  is  known  to  have  a  stability  of  about  110 

kcal/mol  with  respect  to  dissociation  into  HCCN  +  H.^^  The  fraction  of  total  available 
energy  released  as  product  translational  energy  is  0.26  at  both  E^,.  The  remaining  energy 
will  be  partitioned  among  the  vibrational,  rotational  as  well  as  electronic  degrees  of 
freedom  of  the  HCCN  radical  (the  singlet-triplet  splitting  of  cyanomethylene  is  calculated 
to  be  1 1.7  kca/mol).'’”  The  negligible  polarization  of  the  c  m.  angular  distribution  reflects 
a  high  product  rotational  excitation. 

Very  recently,  Takayanagi  and  coworkers'”  have  carried  out  high-quality 
electronic  structure  calculations  on  the  geometry  and  energy  of  the  possible  reaction 
intermediates  and  products  of  the  N(^D)+C2H2  reaction.  Their  work  corroborates  our 
finding  and  in  particular  shows  that  N(^D)  attacks  the  two  C-atoms  of  acetylene  forming 
a  cyclic  intermediate  which  can  lead  to  the  very  stable  CH2CN  intermediate  radical 
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complex  by  (1,2)  H-shift  and  ring  opening,  or  to  a  somewhat  less  stable  linear  HNCCH 
intermediate  by  simple  ring  opening.  The  internally  hot  complex  will  then  dissociate 
under  the  single  collision  conditions  of  our  experiment  to  HCCN  +  H  by  C-H  or  C-N 
bond  cleavage. 

In  conclusion,  the  reaction  of  N(2d)  with  C2H2  produces  HCCN  radicals  which 
are  highly  internally  excited  and,  therefore,  that  may  be  optically  active. 

A  detailed  publication  is  in  preparation.^* 


N(2d)  +  C2H4 

The  reaction  of  N('^S)  with  is  known  to  be  very  slow  (K298=10'^^  cc 

molec‘1  s'^).'*^  Using  beams  of  N(4s,^),  we  have  investigated  the  reaction  with  C2H4 
at  a  collision  energy  of  7.95  kcal/mol  and  detected  products  only  at  m/e=40  (see  Fig. 
21).  This  is  believed  to  originate  from  the  reaction  N(2d)+C2H4  and  to  correspond  to 
the  H-displacement  channel  leading  to  the  formation  of  CH3CN  (acetonitrile)  (AH°o=  - 
107  kcal/mol)  (or  CH2CHN),  which  fragments  completely  in  the  ionizer  to  give 
CH2CN'''  (m/e=40).  Data  analysis  in  terms  of  an  uncoupled  center-of-mass  angular 
distribution  and  translational  energy  distribution  gives  an  average  fraction  of  total 
available  energy  released  in  translation  of  0.23,  which  indicates  that  the  cyano-product 
being  formed  is  highly  internally  excited.  The  c.m.  angular  distribution  shows  a  broad 
forward  peaking  which  suggests  the  intermediacy  of  an  osculating  complex.  Noteworthy, 
acetonitrile  has  been  recently  observed  in  the  atmosphere  of  Titan,  where  both  N(^D)  and 
G2H4  are  present."*^ 

The  experimental  data  could  also  be  fit  assuming  that  molecular  hydrogen 
elimination,  leading  to  CH2-CN  (cyanomethyl)  formation  (AH°p=  -  89  kcal/mol),  takes 
place.  However,  the  resulting  fraction  of  energy  channeled  in  translation  would  only  be 
13%,  a  rather  modest  (and  perhaps  unlikely)  fraction.  We  think  that  formation  of 
CH2CN+H2  is  unlikely,  since  the  three-four  center  elimination  process  which  would  be 
involved  is  expected  to  have  a  high  energy  barrier.  Electronic  structure  calculations  on 
the  energetic  of  all  possible  intermediates  and  reaction  products,  as  those  recently  carried 
out  for  N(^D)+C2H2,'*^  would  be  very  valuable  to  assist  the  interpretation  of  the  crossed 
beam  data.  These  calculations  may  become  available  in  the  near  future. 


Second  RF  discharge  source 

The  construction  of  a  second  radio-frequency  discharge  beam  source,  similar  in 
characteristics  to  that  described  in  ref  20  and  used  in  all  the  studies  reported  here,  has 
been  completed.  The  construction  of  a  new  set  of  RF  electronics  is  at  the  final  stage. 
This  will  allow  us  during  the  next  future  to  attempt,  by  using  the  two  RF  discharges 
simultaneously,  the  dynamical  study  of  radical-radical  reactions  as,  for  instance,  N  +  OH 
->NO  +  H. 
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CONCLUSIONS 

From  product  angular  and  velocity  distribution  measurements  in  crossed  beam 
experiments  we  have  elucidated  the  dynamics  of  the  following  reactions : 

1)  OH  +  CO  ^  CO2  +  H 

2)  OH  +  H2  H2O  +  H 

3)  0(^D)  +  H2  ->  OH  +  H 

4)  0(¥/D)  +  CH3l  ^I0  +  CH3 

5)  N(^D)  +  H2^NH  +  H 

6)  N(^D)  +  C2H2  HCCN  +  H 

7)  N(^D)  +  C2H4  ^  CH3CN  +  H 

The  experimental  results  have  been  compared,  whenever  possible,  with  state-of-the-art 
dynamical  calculations  (quasiclassical  and/or  quantum)  on  ab  initio  potential  energy 
surfaces. 

Reactions  (l)-(3),  and  presumably  also  reactions  (5)  and  (6),  occur  in  the  local 
atmosphere  of  the  orbiting  Space  Shuttle.  It  has  been  shown  that  the  products  of  some  of 
these  reactions,  as  the  OH,  NH,  and  HCCN  radicals,  and  the  H2O  molecule  are  formed 
highly  vibrationally  excited  and  therefore  they  will  give  rise  to  strong  emission  in  the 
infrared.  Furthermore,  the  subsequent  reactivity  of  these  reaction  products  will  be 
considerably  affected  by  their  high  internal  energy  content. 

The  information  obtained  from  the  study  of  these  elementary  reactions  at  the 
microscopic  level,  besides  the  interest  from  a  fundamental  point  of  view,  is  particularly 
relevant  for  a  detailed  understanding  and  modeling  of  the  local  atmosphere  of  the 
orbiting  Space  Shuttle.  In  particular,  it  may  be  useful  for  the  improvement  of  the 
SOCRATES  code,  aimed  at  assessing  the  effects  of  contamination  on  measurements 
aboard  spacecrafts  in  low  Earth  orbit.  Ultimately,  it  should  enable  the  world  scientific 
community  and  the  U.S.  Air  Force  to  better  understand  the  local  environment  of  space 
vehicles  and  of  measuring  platforms  in  space. 
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Figure  Captions 

Fig.  1.  OH  (OD)  product  laboratory  angular  distributions  for  the  0(^D)+H2  reaction 
at  Ec=0.88,  1.9  and  3.0  kcal/mol  (top  three  panels)  and  0(^D)+D2  at  Ec=5.3 
kcal/mol  (bottom  panel).  Solid  lines  ;  calculated  curves  with  best  fit  center-of- 
mass  angular  and  translational  energy  distributions.  Dashed  lines  :  QCT 
predictions  on  ground  state  DIM  surface.  Canonical  Newton  diagrams  are  also 
shown  on  the  r.h.s.;  there  the  circles  delimit  the  maximum  speed  that  the 
OH(OD)  product  can  attain  assuming  that  all  the  available  energy  is  channeled 
into  translation. 

Fig.  2.  Center-of-mass  product  angular  distributions  for  the  same  reactions  of  Fig.  1. 
An  increasing  backward  bias  can  be  noted  with  increasing  Ec. 

Fig.  3.  Schematic  energy  level  and  correlation  diagram  for  0(^D)+H2. 

Fig.  4.  Comparison  between  the  experimentally  derived  OD  product  angular 

distribution  (solid  line)  from  the  0(*D)+D2  reaction  at  Ec=5.3  kcal/mol  and 
QCT  calculations  (dashed  lines)  on  the  ground  state  I’A’  and  first  excited 
state  lA”  PES.  The  summed  contribution  (1  *A’  +  lA”)  is  also  shown. 

Fig.  5.  The  singlet  (solid  line)  and  the  triplet  (dashed  line)  potential  energy  surfaces 
(schematic)  for  the  reaction  0(^P/D)  +  CH3I  ^  10  +  CH3  with  intersystem 
crossing  (ISC). 

Fig.  6.  Laboratory  angular  distribution  (dots)  of  the  10  product  at  Ec=15.3  kcal/mol 
from  the  reactions  0(^P/D)  +  CH3I  ^  10  +  CH3  The  circles  in  the  velocity 
vector  (“Ne\vton”)  diagram  delimit  the  maximum  velocity  that  the  10  product 
can  attain  on  the  basis  of  energy  conservation  if  all  the  available  energy  for  the 
triplet  and  the  singlet  reactions  goes  into  product  translational  energy.  The 
separate  contributions  to  the  total  lab  angular  distribution  (heavy  solid  line) 
from  the  singlet  reaction  (light  solid  line),  the  direct  triplet  reaction  (short 
dashed  line)  and  the  triplet  reaction  via  ISC  (long  dashed  line)  are  shown. 

Fig.  7.  As  in  Fig.  6,  but  for  Ec=13.2  kcal/mol. 

Fig.  8.  Time-of-flight  (TOF)  distributions  of  the  10  product  at  Ec=13.2  kcal/mol.  The 
separate  contributions  to  the  TOF  spectra  are  shown  (notation  as  in  Fig.  6). 

Fig.  9.  Center-of-mass  product  angular  and  translational  energy  distributions  at 
Ec=13.2  kcal/mol. 

Fig.  10.  Beam  transmittance  for  the  supersonic  beam  of  *'’N  seeded  in  Ne  as  a 

function  of  the  magnetic  field  B  at  the  indicated  velocity.  Dotted-dashed, 
dashed  and  long-dashed  lines  are  calculated  assuming  the  exclusive  population 
of '‘S,  ^D  and  ^P,  respectively  ;  the  best-fit  calculation  is  reported  as  solid  line 
and  corresponds  to  a  total  concentration  of  the  two  metastable  states  equal  to 
28%  (^S  :  72%,  ^D  :  21%,  ^P  :  7%). 

Fig.  1 1 .  Beam  transmittance  I/Io  for  the  O2  beam  seeded  in  Ne,  measured  as  a  function 
of  magnetic  field  B  at  v=1210  m/s.  The  long-dashed  line  is  the  calculated  I/Io 
for  a  pure  state,  while  the  dashed  line  refers  to  a  pure  *Ag  state,  both 
calculations  being  performed  assuming  a  rotational  temperature  Tr=600  K  ;  the 
solid  line  is  the  best  fit  calculation  which  gives  the  lower  limit  (-15%)  of  ^Ag 
concentration.  Similar  calculations  assuming  Tr=30  K  give  the  upper  limit 
-40%. 
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Fig.  12  Lab  angular  distribution  of  the  ND  product  from  the  N(^D)+D2  reaction  at 
Ec=3.8  and  5.1  kcal/mol.  Solid  lines  :  calculations  with  best  fit  c.m.  angular 
and  translational  energy  distributions. 

Fig.  13  TOF  distributions  of  the  ND  product  at  selected  lab  angles  for  the  two 
collision  energies  of  Fig.  12. 

Fig.  14  Best-fit  c.m.  ND  angular  distribution  at  Ec=3.8  kcal/mol  compared  with  the 
results  of  QCT  calculations  by  Schatz  et  al.  on  a  new  ab  initio  ground  state 
NHz  PES. 

Fig.  15  Schematic  (semi-quantitative)  energy  level  and  correlation  diagram  for  the 
N(^D)+H2  reaction  based  on  ab  initio  calculations  of  the  ground  and  excited 
PES  by  Harding,  Schatz,  Rabitz  and  coworkers. 

Fig.  16  Energy  levels  for  the  0(^D)+D2  and  N(^D)+D2  reactions. 

Fig.  17  Lab  angular  distribution  of  the  OD  product  from  the  0(’D)+D2  reaction  at 
Ec=5.3  kcal/mol  (l.h.s.)  compared  to  that  of  the  ND  product  from  the 
N(^D)+D2  reaction  at  Eo=5.1  kcal/mol  (r.h.s.).  The  corresponding  Newton 
diagrams  are  also  shown  underneath  the  angular  distributions. 

Fig.  18  Schematic  energy  level  and  correlation  diagram  for  the  N('*S,  ^D,  ^P)  +  C2H2 
system. 

Fig.  19  Laboratory  angular  distribution  of  the  HCCN  product  (detected  at  m/e=38) 
from  the  N(^D)+C2H2  reaction  at  Ec=9.5  kcal/mol.  The  circle  in  the  canonical 
Newton  diagram  delimits  the  maximum  speed  that  the  HCCN  product  can 
attain  if  all  the  available  energy  is  channeled  into  translation. 

Fig.  20  TOF  distributions  at  selected  lab  angles  of  the  HCCN  product  from  the 
N(^D)+C2H2  reaction  at  Ec=9.5  kcal/mol. 

Fig.  21  Lab  angular  distribution  and  TOF  spectra  at  selected  angles  of  the  CH3CN 
product  (detected  at  m/e=40)  from  the  N(^D)+C2H4  reaction  at  Ec=7.95 
kcal/mol. 
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